Dynamical Mean-Field Theory of Resonating Valence Bond Antiferromagnets 
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We propose a theory of the spin dynamics of frustrated quantum antiferromagnets, which is 
based on an effective action for a plaquette embedded in a self-consistent bath. This approach, 
supplemented by a low-energy projection, is applied to the kagome antiferromagnet. We find that 
a spin-liquid regime extends to very low energy, in which local correlation functions have a slow 
decay in time, well described by a power law behaviour and ui/T scaling of the response function: 

x"H «^7(u/r). 



The possibility of frustrated quantum antiferromag- 
nets (QAF) having a resonating valence bond (RVB) 
ground state, that is, a superposition of states where all 
spins are paired into singlets, was suggested many years 
ago such a ground state has neither long-range spin 
order (eg Nccl order) nor spin-Pcicrls order (an ordered 
arrangement of these singlet pairs). Recently, convinc- 
ing evidence has been given that some frustrated two- 
dimensional QAF indeed have RVB physics . In par- 
ticular, numerical studies 0] of the spin-half Heisenberg 
QAF on the kagome lattice reveal that the model has no 
long-range order and displays a very small gap to mag- 
netic (triplet) excitations (estimated to be of order J/20). 
Moreover, this gap is filled with an exponential number 
of singlet excitations, suggesting a possible continuum in 
the thermodynamic limit, and a large low temperature 
entropy, in agreement with the RVB picture. We also 
note, though this is seldom emphasized, that the results 
of P] suggest a correspondingly large number of triplet 
states immediately above the gap. This may be expected, 
since in any valence bond component of an RVB state, 
one can replace any singlet pair by a triplet pair and con- 
tinue to have an eigenstate of S. Resonance between such 
states would lower the excitation energy to much below 
J. 

These considerations suggest that, in a temperature 
range above the triplet gap, the spin correlations have 
a rapid decay in space, but a slow decay in time due to 
the large density of triplet excited states. Equivalcntly, 
the dynamical susceptibility x(q, w) would not display a 
narrow peak around a specific ordering wave- vector, but 
x"{q,w) would have weight at low frequency, character- 
istic of a spin liquid. Indeed, recent inelastic neutron 
scattering studies Q on the S = 3/2 kagome slab com- 
pound SCGO reveal such behaviour above the freezing 
temperature (T > 4K). The only relevant energy scale 
in this spin-liquid regime is apparently set by the tem- 
perature itself, with X\ oc = J2qX"(q, w ) obeying a scaling 
behaviour X\ oc ~ w~ a F{uj/T) with a ~ 0.4, correspond- 
ing to a slow decay in time of the local dynamical cor- 
relations (S(x,Q)S(x,t)) ~ l/t 1 -" ~ l/t 06 . We observe 
that the uniform (<f = 0) susceptibility in this regime is 



well fitted by x oc T~ 0A too. 

In this article, we introduce a novel theoretical ap- 
proach to the spin dynamics in the spin-liquid regime 
of frustrated QAF, which focuses on short-range spin 
correlations only. This approach draws inspiration from 
the dynamical mean-field theory (DMFT) of itinerant 
fermion models j|, and some of its extensions [^|4l0|. 
Our method is fairly general but is applied here to the 
concrete case of the spin-half kagome QAF. The deter- 
mination of on-site and nearest-neighbor dynamical spin 
correlations is mapped onto the solution of a model of 
quantum spins on a triangular plaquette coupled to a self- 
consistent bath. (A single-site approach is not adequate, 
since the essential physics of singlet formation involves at 
least nearest-neighbor sites). Using a projection onto the 
low-energy sector of this model, we demonstrate that, in 
a temperature range extending down to T < J, a slow 
(power-law) decay of temporal correlations is found. 





FIG. 1. The kagome lattice, viewed as an up-pointing clus- 
ter (bold) embedded in a network of similar clusters, which 
is approximated in the DMFT approach by a self-consistent 
bath. 

We view the kagome lattice as a triangular superlattice 
of up-pointing triangular plaquettes (Fig. [l]). Sites are la- 
beled by an index a = 1,2,3 within a plaquette and by a 
triangular superlattice index / numbering the plaquette. 
We denote by Xab(q, t — t') the Fourier transform of the 
dynamical spin correlation function ^{S a j(T) ■ Sb,i'(r')) 
with respect to the plaquette coordinates I, I' (q is a vec- 
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tor of the supercell Brillouin zone). Our approach relies 
on approximating the correlation function by 

X(Q, iv n)ab = J ab{q) + M a b{w n ). (1) 

Here, v n = 2mt/ 1 f3 is a bosonic Matsubara frequency, 
and all quantities are matrices in the internal plaquette 
indices (a, b). J a b(q) is the supercell Fourier transform of 
the exchange couplings. M. a b is a measure of how much 
the correlation function differs from that of a Gaussian 
model, for which \ — J 1 ^ and hence plays the role of 
a spin "self-energy" matrix for the QAF. Obviously, the 
key assumption made in (^) is that the (/"-dependence of 
this self-energy can be neglected. This approximation, 
which is likely to be reasonable when spatial correlations 
are short-ranged, is analogous to the assumption of a 
momentum-independent self energy made in the DMFT 
of correlated fermion models. Here however, the DMFT 
concept is extended on two accounts: the local ansatz 
is made on the spin correlation function rather than on 
a single-fermion quantity, and a plaquette rather than a 
single site is considered. Extensions of DMFT to clusters 
and to (spin or charge) response functions | [n]] 
have been considered separately before in different con- 
texts. Combining them is a unique aspect of our ap- 
proach, which is necessary to capture the dynamical as- 
pects of inter-site singlet formation at the heart of spin- 
liquid behaviour. 

In order to calculate the dynamical self-energy matrix 
■M-ab{iv n ), an effective action involving only the spins of 
a single triangular plaquette is introduced, which reads 

S = S B + ~ J / dr S a - S b + 

i r p 

+ - / dTdT'J2 D °b(r-T')S a (T)-S b (T'), (2) 
Z J ° ab 

in which Sb denote Berry phase terms. D ab (r — r') is 
a retarded interaction, generated by integrating out all 
spins outside the plaquette. Higher order interactions are 
also induced in this process, which have been neglected 
in (||). Equivalently, one can view the rest of the lat- 
tice as an external bath which couples to the spins in 
the cluster through time-dependent external fields with 
a Gaussian correlator D a b- The latter will be determined 
by a self consistency requirement, which stipulates that 
the correlation functions on a plaquette calculated from 



the above action (xab(r - t') = \{Sa{r) ■ S b (r'))s) in- 
volve the same self-energy as the entire lattice. This 
reads: Xab fan) = J&ab + D ab {w n ) + M a b{iv n ), where 
JA a b = J(l — 5 a b) is the matrix of nearest neighbor 
couplings on the cluster. Inserting this relation for the 
self-energy into (Q) and imposing that Y,<?Xab(q) = Xab 
leads to the final form of the self-consistency requirement 

Xab{iv n ) = ^2 [ J( ^ ~~ JA + X~ 1 (wn) - D(iv n )] a l . 

(3) 

We note that the matrix J(q) a b — J^ab involves only 
the exchange constants linking together different (up- 
pointing) triangular plaquettes. In the limit where these 
inter-plaquette couplings vanish while the internal ones 
are kept fixed (decoupled triangles), our approach be- 
comes exact since (||) implies D a b = and (||) reduces 
to the action associated with the Heisenberg model on 
a single triangle. We also note that the above DMFT 
equations can be derived from a Baym-Kadanoff formal- 
ism in which a functional of the correlation function and 
self-energy matrix is introduced in the form 

«[X, M] = ±Trln[M + J] - ^Tv[ x M] + <%]. (4) 

The stationarity conditions |^ = jjj = lead to the 
above equations when the exact functional $ is approxi- 
mated by its value for a single triangular plaquette. The 
free energy of the model can thus be calculated in the 
DMFT approach by inserting the self-consistent values 
of x and A4 into (^) . Alternatively, a functional of the 
correlation function only can be used, along the lines of 
Ref. §. 

In a phase with unbroken translational invariance, the 
matrix Xab{iv n ) in fact reduces to two elements: the lo- 
cal susceptibility xioc = Xaa and the nearest-neighbour 
one Xnn = Xab (a i= b) (and similarly for D a b). It is 
actually convenient to use the following linear combina- 
tions (proportional to the eigenvalues of the \ab and D a b 
matrices): \s = 3(xioc + 2xnn) , Xm = 3 (xioc - Xnn), 
D s = ±(Aoc + 2Am) , D m = | (Aoc - D nn ). In- 
troducing the corresponding self-energies M. s (iv n ) = 
3/xs(^n) ~ 3D s (ii/ n ) , M m (iv n ) = 3/x m (w„) - 
3D m (iis n ), straightforward but tedious algebra allows us 
to recast (0) in the form of two scalar equations: 



X« = f M m + J- § Je 

3 J M s (M m + J) - 2J 2 - \{M. - M m )Je P[e> £ ' 

Xm = f MsM m - J 2 + \J{2M m -M s - J)e 

3 J {M m -J)(M s {M m + J)-2P-l{M s -M m )Jz) 
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Here, e stands for cos q\ + cos q2 + cos #3, where the qi are 
components of q along the three directions of the triangu- 
lar lattice (qi +<Z2+S3 = 0), and p(e) is the corresponding 
density of states. The (^-summation in ([|) simplifies into 
an e- integration in (|^). 

Equations d|J|) entirely define the plaquette-DMFT 
approach introduced here and one could embark at this 
stage into a numerical determination of the two key dy- 
namical quantities XIoc(t) > Xnn(T) using, e.g., a general- 
isation of the quantum Monte Carlo algorithm recently 
introduced in jllj in the context of quantum spin-glasses. 
Instead, we shall gain further insight into the problem 
by making use of a projection onto the low-energy sec- 
tor of the Hilbert space |T^]. Interesting insights into 
the low-energy excited states of the singlet sector have 
been recently obtained by Mila and Mambrini starting 
from the limit of decoupled up-pointing triangles JT^j-lljJ . 
Each such triangle has a four-fold degenerate ground- 
state with spin 5=1/2 and a four-fold degenerate ex- 
cited state with spin S = 3/2. Hence the lattice of de- 
coupled triangles has an exponentially large number of 
ground-states, equal to 4 A ' S / 3 . The physical picture of 
|l3| , fL4| is that the degenerate ground state broadens into 
a band as the inter-triangle coupling is turned on. In the 
following, we shall project the DMFT equations onto the 
low-energy Hilbert space corresponding to the S = 1/2 
sector of each triangle. This amounts to first diagonal- 
izing the single-plaquette hamiltonian (neglecting D a t) 
and then rewriting the retarded term in the action (Q) in 
this low-energy sector. To this end, it is convenient to 
follow Jl2|,[l3| and choose as a basis of the 4-fold degen- 
erate ground-state the eigenvectors of the two following 
operators: the spin- 1/2 operator corresponding to the to- 
tal spin S in a plaquette, and a doublet of Pauli matrices 
T = (T x ,T y ) with the meaning of a plaquette chiral- 
ity operator. In terms of those, the low-energy projec- 
tion of the three spin operators in a plaquette read (with 
lu = expi27r/3 and for a = 1,2,3) 

S 1 ™ = -S [1 - 2u 1 - a T~ - 2uj a - 1 T + ] (6) 
3 

The low-energy projection of the action (||) is then easily 
obtained as 

Slow = S B + drdr' [-D s {t - r') + 
Jo * 

+D m (r-r')f(r)-f(r')} S(r) • S(r% (7) 

where D s and D m are as defined earlier. It turns out 
that Xs and Xm have simple interpretations as the "total 
spin" and "mixed" correlation functions: 

Xs(t-t') = ±{S(t)-S(t')) 

Xm(r-r') = ^{S(r)-S(T')f(T)-f(r')) (8) 



To solve the local quantum problem defined by (Q), 
we use an approximate technique that has proven suc- 
cessful in recent studies of quantum spin-glass mod- 
els pp[ and of impurities in quantum antiferromagnets 
p5[ . We introduce two doublets of spin- 1/2 fermions 
(sj,sjj and (if ,£,[.) subject to the constraints: s^s-f + 
s^si = t^t-^ + tfti = 1 in order to represent the op- 
erators S and T as: S + = s^sj. , S~ = s+Sj , S z — 
(s+s T - s+s^/2 - T+ = t+t i ,T~ = t^tt. The interact- 
ing fermion problem corresponding to ([?]) is then solved 
in the self-consistent Hartree approximation. Also, the 
local constraint is approximated by its average (the as- 
sociated Lagrange multiplier is zero due to particle-hole 
symmetry). Introducing the Green's functions G s and 
Gt for the two fermion fields and defining self-energies by 
G7 x (w n ) = iuj n - E s (iu n ) , G^ x (iw n ) = iuj n - S t (iw n ) 
(with uj n = (2n + l)ir/f3 a fermionic Matsubara fre- 
quency), this yields the imaginary-time equations: 

S s (r) = - ^D s (t)G s (t) + 3D m (r)G s (r)G t (r)G t (-r) , 
Et(r) = 3D m (T)G t (T)G s (T)G s (-T). (9) 

Then the local problem (for a given bath (D s ,D m )) re- 
duces to two coupled non-linear integral equations. The 
correlation functions are given by 

Xs{t) = -\g s {t)G s (-t), 

Xm{r) = 2G s (r)G s (-T)G t (r)G t (-r). (10) 

In practice, we use the following algorithm: we start 
with an initial guess for the bath D Si m, and obtain the 
Green's functions by iteration of e q. (M ). We can then 
calculate the susceptibilities from (|lC|), as well as the 
self energies M s , m - Inserting the latter into (||) yields 
new Xs,m and, in turn, the new baths D^™(iiy n ) — 
-M St m(iv n )/3 + t/Xs,m{iv n )- 

We now discuss our findings when solving these cou- 
pled equations numerically. The first key observation is 
that a paramagnetic solution can be stabilized down to 
the lowest temperature we could reach, with no sign of 
long-range ordering intervening. Long-range order is as- 
sociated with a diverging eigenvalue of Xab(q) for some 
q and hence, within our scheme, to a vanishing denomi- 
nator in ([5]), which we never observe. In Fig. ^, we dis- 
play our results for the temperature dependence of the 
local (i.e on-site, or q- integrated) susceptibility xioc{T) 
— L x\oc(T)dr. At high temperatures xioc oc 1/T obeys 
Curie's law. Below T « 0.5 J, the effective Curie constant 
decreases with decreasing T, indicating gradual quench- 
ing of the local moment. However, xioc itself continues 
to increase as the temperature is lowered. In the inset 
of Fig. ^, we display the effective exponent defined by 
a(T) = — lnxi oc /lnT. A transient regime, correspond- 
ing to a slowly decreasing a, is apparent and extends over 
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more than two decades. At the lowest temperatures, a 
appears to saturate to a value a ~ 0.5. 
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FIG. 2. Main plot: xioc obeys Curie law (oc 1/T) at high 
temperature, but diverges more slowly than 1/T at lower tem- 
perature. Inset: Evolution of the effective exponent a from 
high to low temperature when fitting xioc(T) to 1/T a . (T is 
in units of J in all figures.) 




FIG. 3. For 

different (low) temperatures, Xioc(t)/xioc(/?/2) plotted as a 
function of r//3 collapse on a single curve, well described by 
[sin(7rr//3)]-°- 5 (solid curve ) . 

Correspondingly, the correlation function Xioc( T ) 
obeys a power-law dependence on time in the low- 
temperature regime. Figure ||shows xioc( T )/xioc(/3/2) as 
a function of r//3 for several (large) values of (3 = 1/T; 
the graphs collapse on each other for different (3, and 
are well fitted by the expression X1oc(t)/xioc(/3/2) = 
[sin nr/j3] ^ , with 1 — a ~ 0.5. This scaling function 
is the one appearing in the study of quantum spin-glasses 
|io| and quantum impurity models and dictated 

in the latter case by conformal invariance. It corresponds 
to the following scaling form of the dynamical suscepti- 
bility: X £» <x Lu- a F a (iu/T) with F a (x) = x a \T(^ + 
* T~ ) 1 2 Sm ^ § • This also implies a diverging behaviour of 
the relaxation rate 1/Ti oc T X [' oc (lo)/lu ~ l/T a . 

These findings can be interpreted as the formation of 



a spin-liquid regime with a large density of triplet ex- 
cited states, consistent with ref. ||. Our approach fails 
to predict a triplet gap, however, because our triangular 
plaquette has a spin-1/2 ground-state, and the classical 
bath to which it is coupled cannot fully screen this local 
moment. Hence the local susceptibility continues to in- 
crease at low temperatures in our approach, down to an 
unphysically low energy scale (of order J/100). In fact 
the triplet gap for the 5=1/2 kagome QAF is known Q 
to be quite small (of order J/20), so our description of the 
spin dynamics should be valid down to a rather low tem- 
perature scale. The power-law behaviour of the spin cor- 
relations and the low-temperature increase of the suscep- 
tibility agree well with experimental findings on SCGO 
above its freezing temperature at T = 4K g . Since this 
is a S = 3/2 system, it is conceivable that higher spin 
extends the range of applicability of our approach even 
further. A promising application of our approach is the 
pyrochlore lattice: the natural plaquette is a tetrahedron, 
with a twofold degenerate singlet ground state, so that 
the very low-energy singlet sector should be accessible 
within our approach. In future work we plan to address 
these issues, and also study the low-temperature thermo- 
dynamics by solving the self-consistent local problem (||) 
using exact numerical techniques. 
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